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Abstract 

Understanding of the signatures of cosmic string networks is limited by a large 
uncertainty in the sizes at which cosmic string loops form. We review cosmic string 
network evolution, and the gravitational signatures, with emphasis on this uncertainty. 
We then review a recent analytic model of cosmic string networks. In combination with 
recent simulations, this suggests that 90% of the string goes into very small loops, at 
the gravitational radiation scale, and 10% into loops near the Hubble scale. We discuss 
cosmic string signatures in such a scenario, and the 'inverse problem' of determining 
the microscopic cosmic string properties from observations. 

To appear in the Proceedings of the 11th Marcel Grossman meeting. 



1 Introduction 



Compactifications of string theory give rise to many potential cosmic strings, including 
the fundamental strings themselves, D-strings and wrapped D-branes, solitonic strings and 
branes in ten dimensions, and magnetic and electric flux tubes in the four- dimensional ef- 
fective theory. In various models of inflation in string theory or grand unified theories a 
network of such strings will form. Thus it is interesting to ask, what are the prospects for 
future discovery of these strings if they exist, and, what are the prospects for distinguishing 
different microscopic models'!]]] 

In between the microscopic theories and the observations, the strings have a variety of 
macroscopic parameters and properties, including: 

• Their tension, /x. 

• Their intercommutation probability P. When two strings collide, they either pass 
through each other, with probability 1 — P, or intercommute (reconnect), with proba- 
bility P. In fact, P is a function of the velocity and angle of the collision, but to first 
approximation one can consider the value averaged over the kinematics. 

• Their gapless degrees of freedom: are these just the collective coordinates for the motion 
in three-dimensional space, or are there additional bosonic or fermionic modes? 

• Their interactions with low energy fields: are these gravitational only, or axionic or 
gauge as well? 

• The number of kinds of string: one or many? 

• The existence, or not, of multi-string junctions. 

One would hope to proceed from the observation of cosmic strings to a determination of 
their macroscopic properties, and from there to possible microscopic models]^] This requires 
the ability to calculate with some degree of precision the signatures of a given macroscopic 
model. In fact, this is quite challenging. 

Here we will focus on the most vanilla type of cosmic string: a single species with P = 1, 
no long-ranged interactions besides gravity, and no multi-string junctions or extra degrees 
of freedom. Even for these, major uncertainties remain. Fig. 1 shows a simulation of the 

1 My MG11 talk had substantial overlap with hep-th/0410082 [1 , to which the reader is referred for 
further details and references. Here I will concentrate on some recent developments, in particular my work 
with Jorge Rocha and Florian Dubath. This is a slightly longer version than in the official proceedings, 
with the addition of Sec. 3.4.2 and related discussion in Sees. 4.2 and 5.2, and a new title (the title in the 
proceedings, 'The Cosmic String Inverse Problem' seemed too broad in retrospect, since we consider almost 
exclusively vanilla strings) . 

2 That is, we are separating the subject into 'model building,' the connection between the microscopic and 
macroscopic models, and 'phenomenology,' the connection between macroscopic models and observations. 
The present paper (and our references) will be focussed entirely on phenomenology. 
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network, with one side of the box being about a quarter of the Hubble length. This simulation 
is for the radiation era; the matter era simulations have the same general features but some 
quantitative differences. One sees several long strings crossing the simulation volume, and a 




Figure 1: Simulation of a vanilla string network during the radiation era, from Allen and 
Shellard [3] . The scale of the box is approximately one quarter of the Hubble length. 

gas of string loops. 

As will be explained in more detail in the next section, the original expectation was 
that the Hubble length (which is the same order as the horizon length, and the FRW time) 
would be the only relevant scale in the problem [2J. If this were true, then simulations, 
beginning with Refs. [31 ffl [2], should have readily produced good quantitative descriptions 
of the networks. However, inspection of Fig. 1 shows that the long strings are not smooth 
at short distance, but have a kinked structure, which persists down to the cutoff scale of 
the simulations. Further, it appears that the loops are rather smaller than the Hubble scale. 
Indeed, in the early simulations the loop distribution seemed to peak near the short distance 
cutoff scale of the simulations. 

Subsequent work has failed to produce a clear picture, and in fact the estimates of the 
typical size at which the loops are produced range over more than fifty orders of magnitude. 
This is rather remarkable, when one considers that the entire span from the Planck scale 
to size of the visible universe is just sixty orders of magnitude. Moreover, the problem is 
completely well-posed: the limitation is only in our ability to solve the equations. The 
problem with the simulations is that they are limited to just a few orders of magnitude 
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in length scale, and even more limited in time scale: for the latter not only is there the 
direct cost, but the expansion of the universe means that increased time also costs increased 
simulation volumej^] Much of the work is based on an expansion factor that is only of order 
5, and it is difficult to distinguish transient effects from real ones. 

Analytic methods are difficult because of the nonlinearities of the problem. One can 
write down the evolution equations for an ensemble of strings [6], but before long one must 
begin to approximate, and in the end most analytic treatments reduce to somewhat coarse 
models in terms of a few parameters. On the other hand, one might have expected that the 
nature of the problem, the existence of a large ratio of scales, would lend itself to something 
like a renormalization group treatment. Roughly speaking, one would use the simulations to 
treat the strongly nonlinear behavior at the Hubble scale, and then evolve down to shorter 
scales analytically. 

This was the motivation for the recent work of Refs. [TJ El [9j [10] . Thus far, things are not 
as systematic as for the renormalization group. The difficulty is evidently that for the string 
network the flow is not from short distances to long but the other way around. The Hubble 
length, against which all things are scaled, increases more rapidly than the comoving length 
(during the matter and radiation eras), so comoving structures move to smaller effective 
scales over time. In this respect the problem is like turbulence, and has been similarly 
frustrating. 

Nevertheless, it has been possible to proceed by making approximations that are not 
precisely controlled but perhaps well-motivated. In particular, this has allowed us to un- 
derstand why structures appear at scales far below the Hubble length. The simulations are 
still essential, but the analytic model enables us to distinguish real effects from transient 
ones, and to extrapolate beyond the scales that can be simulated. Combined with recent 
simulations that use a trick to get beyond the limitation on expansion times, it may be that 
a consistent and fairly quantitative picture is emerging. 

In Sec. 2 we review the basic features of network evolution, defining the problem. In Sec. 3 
we discuss some aspects of the gravitational wave signatures, focussing on the sensitivity to 
loop size. In Sec. 4 we present the analytic model of the short distance structure. In Sec. 5 
we attempt to compare the model with the simulations, and discuss the current picture. 

2 Network Evolution 

In this section we review the basics of network evolution, focusing on the issues that will be 
relevant later. Aside from some recent work, this is standard material, on which the reader 
can find more detail in the excellent book by Vilenkin and Shellard 
The relevant processes for vanilla strings are: 

3 The horizon length is increasing more rapidly than the comoving length, so the number of degrees of 
freedom that can interact with one another increases with time. 
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1. The formation of an initial string network, such that a finite fraction of the string goes 
into infinite random walks. 

2. The stretching of the string network with the expansion of the universe. 

3. The intercommutation of long strings, producing kinks. 

4. The emission of gravitational radiation from long strings. 

5. The self- intercommutation of long strings, producing loops, which may fragment fur- 
ther into smaller loops. 

6. The decay of the loops via emission of gravitational radiation. 

We will discuss each of these processes in turn, and then the key issue of the scaling property 
of the network. 

2.1 The Kibble Process 

Any gauge theory with a broken £7(1) symmetry has a string soliton solution. Moreover, 
such strings will actually form in any phase transition in which a U(l) symmetry becomes 
broken [2j. At the transition, the Higgs field rolls from the origin to one of the symmetry- 
breaking vacua, and unless the transition takes place infinitely slowly, causality implies that 
the phase will uncorrelated over long distances^] Thus defects - strings - will form. Causality 
implies that over long distances the strings form random walks. A random walk has a finite 
probability not to return to its starting point, so some of the string will be in infinite random 
walks; this conclusion is supported by simulations. These infinite strings are the seeds of the 
later network, because any loops that form in the transition will rapidly decay. 

Subsequently, Sarangi and Tye [T2l [T5] made the remarkable observation that D-strings 
would be produced by D-brane annihilation at the end of brane inflation, since the low 
energy dynamics is again described by U(l) symmetry breaking. In fact, regardless of the 
microscopic structure of the strings, there is be some dual version of the Kibble process, such 
that if the strings can exist only after some phase transition, they will actually form during 
that transition. 

Stability 

We must assume that the processes of breakage and string confinement [T4"t [T5| [T6] are either 
absent or sufficiently suppressed so as to be negligible on cosmological time scales. 

4 Of course the phase is not gauge invariant, but we can compare it using parallel transport, or in Coulomb 
gauge. 
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2.2 Strings in Expanding Spacetime 



Cosmic strings are long compared to their thickness, and so can be treated as idealized one- 
dimensional objects. For vanilla strings, the relevant action is just the Nambu action in the 
FRW metric 

ds 2 = —dt 2 + a(t) 2 dx ■ <ix 

= a(r) 2 (-dr 2 + dx • dx) . (2.1) 

The equation of motion governing the evolution of a cosmic string is 

x + 2-(l -x 2 )x = - (— ^ . (2.2) 
a e V e J 

Here e is given by 

x /2 xl/2 

These equations hold in the transverse gauge, where x ■ x' = 0. Dots and primes refer to 
derivatives relative to the conformal time r and the spatial parameter a along the string, 



respectively. The evolution of the parameter e follows from equation (2.2) 



- = -2-x 2 . (2.4) 
e a 

The equation of motion implies that on scales long compared to the horizon length, which 
is of the order of the FRW time t, the string is frozen in the comoving coordinates x, and so 
just expands with the universe. When a scale comes inside the horizon, the modes begin to 
oscillate; the energy begins to redshift away, and the string straightens. 

From the second derivative terms it follows that signals on the string propagate to the 
right and left with da = ±dr/e. Thus the structure on a short piece of string at a given time 
is a superposition of left- and right-moving segments. In an expanding spacetime the left- 
and right-moving waves interact — they are not free as in flat spacetime. In terms of left- 



and right-moving unit vectors p± = x± -x', the equation of motion (2.2) can be written as 



P± =F "P± = [Pt - (P+ ■ P-) P±] i ( 2 - 5 ) 

d (X 

so that in the limit that we can ignore the time derivative of the scale factor we recover 
independent left- and right-moving waves. 

2.3 Long String Intercommutation 

When two straight strings intercommute, the resulting strings have a kink, which separates 
into left- and right-moving kinks. Under the flat spacetime Nambu equations, the kinks 
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would move indefinitely in the two directions. The curved spacetime equations imply that 
the kinks persist but their angle slowly decreases [5]. The decrease is sufficiently slow that 
the kinkiness is potentially very large: if one considers the total root-mean-square kink 
angle in a length of long string (which is what would be relevant if the kink directions were 
uncorrelated) , one finds that it diverges, so there is the potential for having structure on 
very small scales in the network. 

2.4 Long String Gravitational Radiation 

The oscillations of the long string will emit gravitational radiation, and so are damped 
below a certain length scale. A naive estimate would lead to the same scale that we will find 
below for the loop decays, something of order the dimensionless string tension G\i times the 
FRW time t. However, things are more subtle for long strings. Emission of gravitational 
radiation requires that a right- and left-moving mode meet. There is a nonlinear suppression 
when these have very different wavelengths, because the short-wavelength modes essentially 
perceives a straight string [T7], so that the actual damping scale is 0(t) times a larger power 
of Gfj, [HI E]. We will be more quantitative about this in Sees. 4 and 5. 

2.5 Loop Formation 

When a long string intercommutes with itself, a loop breaks off. The loop may then self- 
intercommute and fragment into smaller loops. This happens rapidly, within one period, 
so we are interested in the size distribution of the ultimate non-self-intersecting loops. Of 
course a loop may collide with a long string (or each other) and reattach; for loops close to 
the FRW scale this is likely, but smaller loops have a high probability never to rejoin. 

2.6 Loop Decay 

Dimensionally, a loop of length I emits gravitational radiation power of order Gfi 2 . Since its 
total energy is of order fj,l, we conclude that the lifetime is 

t(l) = l/TGn , (2.6) 

where the numerical factor T is found to be of order 50 for typical loops. In other words, 
loops smaller than TGfit will decay in a Hubble time. 

2.7 Scaling and Loop Size 

Now let us put the preceding ingredients together. A key concept is scaling, the idea that 
the statistical properties of the network should be constant when measured in units of the 
FRW time t [2j. For example, the mean separation between long strings should grow as 
t. If the network simply expanded with the universe, this separation would grow as the 
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scale factor a(t), which increases more slowly than tj^] Scaling requires that the processes 
discussed above reduce the amount of string, and that they do so as fast as causality allows, 
since the horizon distance is also linear in t. For a network that scales, the energy in strings 
in a box of side t is linear in t, so the density is proportional to t/t 3 = 1/t 2 . Thus it is a fixed 
multiple of the dominant energy in each epoch, which also is proportional to H 2 oc t~ 2 . If 
the energy density did not scale, the strings would eventually come to dominate, unless the 
string tension were very low. 

In fact, it appears that string networks scale under rather general conditions, so that 
the scaling solution is an attractor. If the string density exceeds the attractor value, the 
rate of long-string collisions increases, resulting in more kinks and more loop formation; if it 
is smaller, than the opposite occurs. A simple one-scale model captures this behavior [19J. 
Suppose that there is a single length scale L that governs the network, that is, the long-string 
separation, the correlation of structure along the long string, and the formation of loops. 
A convenient definition is in terms of the energy density in long strings, = /i/L 2 . The 
energy of long strings in a comoving volume would then satisfy 

The first term is from the expansion of the universe: the 1 in parentheses is just the growth 
of the length with the scale factor, while the — 2v 2 is from the redshifting away of energy 
below the horizon scale. The latter follows from the definition (2.3), which implies that the 
total energy in a string is simply \i J daae, together with the equation of motion (2.4). The 
second term is from loop formation, with c a constant and the scale L inserted by dimensional 
analysis. 

The average v 2 plays an important role in the properties of the network. If the universe 
were not expanding, the virial theorem would imply that v 2 = 0.5. The expansion of the 
universe reduces this to around 0.41 in the radiation era and 0.35 in the matter era, according 
to the simulations [20J. 

Inserting = a 3 fi/L 2 , and defining x = L/t and a oc t r , one finds that 

— r = (l + v 2 )r-l + — . 2.8 

xdt y ' 2x y ' 

Since v 2 is less than 0.5, the first term is less than 1 in both the matter and the radiation 
eras, and so x is driven to a fixed point value c/2(l — [1 + v 2 ]r). In fact this model is overly 
crude and has been superseded by more elaborate models with multiple scales and other 
parameters, but it captures the basic feedback that is operating. 

The scaling property means that the initial conditions from the formation of the network, 
in particular the initial density of strings, are washed out in time. Thus there are only three 
processes that enter: expansion, intercommutation, and gravitational radiation. However, 



5 Note that this discussion is relevant only during the periods of matter and radiation domination. We 
will discuss some effects of cosmic acceleration later. 
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there are some important questions that remain open: Does the scaling property apply 
only to the gross features of the network, or to all details including the loop distribution? 
And, does gravitational radiation only play the role of removing the loops from the network, 
or is gravitational smoothing of long strings important as well? Note that if long-string 
gravitational radiation can be ignored, then the problem has no parameters: the processes 
of expansion, intercommutation, and loop removal are purely geometric (the string tension 
[A scales out of the equations of motion, and P has been set to 1). 

As we have noted in the introduction, there is still uncertainty about the small scale 
properties of the network. For example, the estimates of the characteristic size Zf at which 
loops form range over 

• If ~ O.lt. A loop size near the Hubble scale was the initial expectation in this subject, 
as embodied in the one-scale model, and in fact this is also the conclusion of some of 
the most recent simulations [2TJ E2] ■ 

• Zf ~ 10~ 2 to 10 3 t. Other recent simulations [221 EH] point to a loop size that still 
scales as a pure number times the FRW time, but rather smaller than the above. 

• Zf ~ TG^it. These loops still scale, but only with the inclusion of gravitational backre- 
action to smooth the long strings [24] . 

• Zf ~ r(G[i) k>1 t. Again this length scales, but as noted earlier backreaction is less 
efficient than originally thought [ITJ [TH1 [S] • 

• Zf ~ string thickness. This length does not scale, but remains at a fixed physical 
size [251 ESI E7]. 

This has been a contentious subject, and indeed in the references a number of names are 
attached to more than one estimate. In the end this will not be surprising, as we will argue 
that more than one of these estimates is correct. 

3 Network Signatures 

Since vanilla strings have only gravitational long-range interactions, the key parameter is 
Gfi, which determines the typical metric perturbation produced by a long string. In models 
of brane inflation, this is related to the scale of CMB fluctuations ST/T, but in a way that 
depends on the geometry so that one obtains as an interesting range [13] 



This is probably also relevant to more general models in which strings are produced at the 
end of inflation. Possible gravitational signatures include 



KT 12 < G/2 < 10 



6 




s 



• Dark matter. 

• Effect on CMB and galaxy formation. 

• Gravitational lensing. 

• Gravitational radiation. 

3.1 Dark matter 

We can immediately eliminate this. We have noted that the string density is a fixed fraction 
of the dominant energy density in a scaling regime. The ratio is proportional to the string 
tension, in the dimensionless combination G[i\ the simulations give approximately 70G/J, in 
the matter era and 420G/Z in the radiation era. Since the current upper bound on G\i is 
around 2 x 10~ 7 , this is rather less than one. Of course, if the network does not scale and 
if the string tension is just right, then the strings could become dominant at some late time 
- but then they would not redshift like matter. The effective value of the equation of state 
parameter would be between and —1/3, so they are also not the dark energy. 

3.2 CMB and Galaxies 

The CMB perturbations come primarily from the long strings, since their fields extend 
coherently over long distances. Thus they depend on network properties that are reasonably 
well understood]^] The scaling property implies that this perturbation spectrum is scale 
invariant, but the power spectrum does not match observations: there are no acoustic peaks, 
because the perturbations are generated after the modes reenter the horizon rather than (as 
with inflation) before they leave. This translates into an upper bound [28], [29] Gfi <2 x 
1CT 7 . Bounds from non-gaussianities are currently less sensitive than those from the power 
spectrum. The strongest bound, around 3x 1CT 7 comes not from limits on string-like features, 
but from limits on Doppler deviations from a black body spectrum that would be introduced 
by the transverse motion of the strings [30] ■ The power spectrum bound is argued to be 
approaching the cosmic variance limit, but bounds from polarization [281 EH] or pattern 
search will improve one or two orders of magnitude with time. 

If the loops are large enough, they can act as seeds for early galaxy formation and so 
lead to early reioniziation [31]. If the largest estimates of loop size are correct, this leads to 
a stronger bound 3 x 10 -8 , with some uncertainties. 

6 However, Ref. [27] finds densities several times lower, so the bounds are several times weaker. It also 
rather optimistically suggests that a best fit per degree of freedom includes cosmic strings. 
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3.3 Lensing 



Both long strings and loops (if they are large enough) can lens more distant objects. The 
deficit angle at a string is 87cGfi, or around 1 arc-second at the current CMB bound. The 
separation for a given object will be somewhat larger or smaller than this, depending on the 
geometry and the string velocity [32]. Occasional lens candidates have been reported, but 
have thus far always proven to be binary galaxies. To obtain bounds from lensing requires a 
large survey, because a string network lenses only a small fraction of the sky. Optical surveys 
may reach Gfi ~ 10 -8 [33], and radio frequency surveys may eventually be sensitive down to 
Gn ~ 10" 9 [33] . 

For lensing, there is the issue of how straight the string is at short distances. Most 
phenomenology has considered two extremes: a straight string, which simply translates 
a strip of the background and so produces identical double images, or a random walk, 
which produces complex multiple images. Further, with the straight string any additional 
lensed objects would be on a line through the lens axis; with the random walk they will be 
distributed more irregularly. It will be one of our results to resolve this for the real networks. 



3.4 Gravitational Radiation 

This seems the most likely route by which cosmic strings will be discovered, or the bounds 
tightened to an interesting degree. Again, we refer the reader to Vilenkin and Shellard [UJ 
for a detailed treatment, as well as Ref. [35], and also Refs. [361 E3 EH] for recent updates. 
In contrast to the signatures above, the gravitational radiation comes predominantly from 
the loops, because these are smaller and so have higher frequencies. Thus the issue of loop 
size is particularly important. 

For vanilla strings, all of the energy going into loops ends up as gravitational radiation. 
We can get the total energy going into loops from energy conservation. As in Sec. 2.7, we 

f[l + «W-^p. ( 3 - 2 ) 

,^l-[l + v 2 }r) Poo . (3.3) 

There are two key distinctions to make: between the radiation from the low and high har- 
monics on a loop, and between loops that are large or small compared to TGpt. Smaller 
loops decay essentially at once on the cosmological time scale, so this energy goes directly 
into gravitational waves. Loops that form at larger sizes persist for a while before decaying. 
We consider the latter case first. 



nave 

dpoo _ 
dt ~ ' 

Putting in scaling, oc t~ 2 , gives 

^Ploop 

dt 
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3.4.1 Large Loops, a > TGjj, 



Once a loop forms (below the Hubble scale), it stops redshifting, so the energy density in 
loops behaves like matter. During the matter era, this is not so interesting, but during the 
radiation era it implies that the energy density in a set of string loops grows relative to the 
dominant energy p r . A large loop that forms at size oetf will decay at td, where at{ = TGptd- 
Thus there is a relative enhancement by a(td)/a(tf) = (a/TG/i) 1 / 2 , if the loop both forms 
and decays during the radiation era. For these loops, then, 

t/-^ = 2(1 - [1 + v^j-^ia/TG^^ . (3.4) 
aid o p T 

We have included two additional factors found in a more complete treatment: a | from 
integrating the a(td)/a(t) over the life of the loop, and a factor of 7~ 3 / 2 where 7 is the 
Lorentz factor with which the loop is produced. A factor of 7 _1 enters directly because the 
kinetic energy redshifts away rather than becoming gravitational radiation, and an additional 
7~ 1//2 enters through a{td)/a(t{) because this reduces the lifetime of the loop. 

A loop of length / has period 1/2 and fundamental frequency v = 2/1. Most of the energy 
goes into the lowest few harmonics, so v{td) ~ 2/TGfitd- The observed frequency is 

a(t d ) ( . a(t d ) 2 

The condition that this be during the radiation era is Gfi > 1.5 x 10~ 10 / (uq • yr). For LIGO 
(z/ ~ 10 2 Hz) and LISA (z/ ~ 10~ 2 Hz) this inequality is satisfied down to very small 
tensions. For pulsar timing (z/ ~ 0.1/yr) it is satisfied down to G\i around 10~ 9 (for smaller 
values of G/j,, the gravitational radiation in the pulsar range first rises and then falls steeply). 

Finally, noting that during the radiation era aV /z/o = —dtd/2td, and that the gravita- 
tional wave energy redshifts like radiation, we have 

^ = 4(1 - [1 + v*]r T ) \r i/2 (a/TG^ 2p -^n T 

= 4(1 - [1 + v 2 r }r T ) L-y^a/TGp) 1 / 2 ^! + z eq )^fi m 

= 0.0035 7 " 3/2 (aG/i) 1/2 . (3.6) 

In the last line we have inserted v 2 = 0.41, r r = ^, Poo/p r = 420G/1, z eq = 3200, and 
Q m = 0.24. Note the factor of Q m : the fractional energy density in gravitational waves is 
diluted by the vacuum energy: this is a suppression relative to older estimates. 

The observed regularity of pulsar signals puts limits on the spacetime fluctuations through 
which they pass, such that [39j i/odficw/^o < 4 x 10~ 8 for z/ ~ (10 yr)^ 1 . Thus 

Gp < 1.3 x 10- 10 crV . (3.7) 
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If we take the largest estimate, a ~ 0.1, and assume that j m 1, we get the strong limit 
G/i < 1.3 x 10~ 9 . The pulsar bounds are likely to improve substantially with time. Further, 
Advanced LIGO is expected to reach values of i/o^cw/^o between 10~ 8 and 10~ 9 , and 
LISA will reach 10~ 12 or better. This makes clear the importance of understanding the 
loops. 

The emission of gravitational radiation from the high harmonics of the loops is larger 
than it would be if the string were smooth. This is first of all due to the kinks, which give 
a z/~ 5//3 spectrum (z/~ 2 / 3 per log frequency), but an even more interesting effect is due to 
the cusps. The functions p± represent two curves on the unit sphere, which will generically 
intersect at points [30] . Each intersection represents a periodic event, in with the string snaps 
like a whip whose tip approaches the speed of light. This produces a cone of gravitational 
waves, with a power spectrum i/ -4 / 3 (z/ -1 / 3 per log frequency) and a distinctive signature in 
time. 

Initial calculations of this signal [UJ H2] indicated a good chance that it would be seen 
at LIGO, but more recent calculations [15] . which treat the cosmology and the detector 
characteristics more precisely, indicate that we may have to wait for LISA. The cusp signal 
is less sensitive to the loop size (as long as it is at least TG/xt) because over much of parameter 
space one is seeing cusps from the matter era. However, the largest loop size, I ~ O.lt, gives 
a signal within the reach of Advanced LIGO down to Gil ~ 10~ 9 [15] . 



3.4.2 Small Loops 

It is customary to write the formation sizes of these in terms of a parameter e < 1, where 
If = eTGiit. For these, the decay is essentially instantaneous and so the rate of gravitational 
wave production is 

t d ^ = 2(l-[l + v 2 m }r m ) Poo . (3.8) 
For the low frequency loops, we have 

* . (3.9) 
ao el Gyutd 

Observations at a given frequency v$ are therefore sensitive only to the parameter range 

Noting that f^cw will be diluted by a factor of f2 m a(id)/a.oQ we have for the low-harmonic 
radiation 

z/ — = 6(1 - [1 + v m \r m ) it m . (3.11) 

au a p m 



7 For loops that form during the accelerating era one needs to redo the network models, but dilution by 
f2 m should be a good estimate. 
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On the other hand, for high-harmonic radiation (log-frequency spectrum v 1//3 ) emitted in 
the very recent past we have 



GW 



dun 



6c(l- [l + v 2 m ]r m )(2/eTGfiu to 



,V3 P^Q r 
Pm 



(3.12) 



Here c is a possible suppression factor if many loops do not have cusps. In Ref. [9] it is 
argued that small loops are likely to have cusps, so this factor is probably close to 1. In 
the range (3.10), for c of order 1, the high-harmonic signal (3.12) is greater than the low- 
harmonic signal (3.11 ): the suppression of the latter due to redshifting is significantly greater 
than the suppression of the former from the falloff of the spectrum (and this would still be 
true even if the high-frequency spectrum were only due to kinks). Inserting the appropriate 
constants, Eq. (3.12) becomes 



^0—7— = 3.5c(ez/ t ) 



1/3 



du t 



(Gfi) 



2/3 



(3.13) 







The bounds (3.10,3.13) are in approximate agreement with the sensitivity found in Ref. 
(their Fig. 2, extrapolated to P = 1). 



Note that Eq. (3.13) represents the contribution of the cusps when they are averaged out 
into a stochastic signal. In Sec. 5.2 we will discuss the observation of individual cusps. 



4 A Model of Short Distance Structure 
4.1 The Long-String Two-Point Function 

We now attempt to get an analytic understanding of the small-scale structure in the net- 
work |8]. Thus we focus on the evolution of a short left- or right-moving segment on a 
long string. This will potentially involve the following: 



1. Evolution according to the FRW-Nambu equations Eqs. (2.5). 

2. Long string intercommutation. 

3. Incorporation into a loop larger than the segment. 

4. Emission of a loop comparable to or smaller than the segment. 

5. Smoothing via gravitational radiation. 

The probability of the second of these is proportional to the length of the segment, and 
so can be systematically neglected for a short segment. The third process is governed by 
dynamics on a longer distance scale, and so does not depend directly on the configuration 
of the segment: it will not change the statistical distribution for configurations of small 
segments. 
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In saying this we are neglecting the correlation between the structure at short distance 
and the structure at long distance. This is not a fully controlled approximation, but should 
become progressively better at shorter scales. One could attempt to make an improved model 
in which this correlation is parametrized, but this will not be necessary for our purposes. 
Ultimately one would like to write down the exact equations for an ensemble of string 
networks, and then approximate systematically, but this is apt to be very difficult; it is more 
practical to first figure out what is going on. 

The small loop production would be absent if the only scale were the Hubble length. 
We will see that if we assume this and proceed, it will not be self-consistent. However, we 
will also find that the production of small loops is primarily controlled by the long-distance 
structure rather than the short distance structure, so we seem to get the right answer anyway. 
Finally, we will assume that, while we are below the Hubble scale, we are at large enough 
scales to ignore gravitational radiation. 

Thus we need only to solve for the motion of a string in an expanding universe. We recall 
the equations for the unit vectors p± = x ± *x', 

P± T -V± = ~ [Pt - (P+ ■ P-) P±] • (4- 1 ) 

These equations are nonlinear and not solvable in closed form, but again we can simplify for 
a short segment. We separate the configuration into a mean direction of the segment and a 
fluctuation, 

P±(r,a) = P ± (r)+w ± (r,a)-^P ± (r)«;|(r,(T) + ... , (4.2) 

where P± = 1 and P± ■ w± = 0. We will expand in powers of w±(r, a): the deviation from 
the mean decreases with the size of the segment (one can show this from the full equations, 
and we will confirm that it is self-consistent). Then one finds 

P + = --[P_-(P + -P_)P + ] , (4.3) 

(X 

W+ - -w' = -(w+ ■ P+) P+ + -CP+ ■ P_) w+ . (4.4) 
e a 



In Eq. (4.4) we have dropped a term ^(P+ • w_)P + , because it averages to zero (with 
corrections suppressed by the length of the segment) as the left- and right-moving waves 
sweep past each other. The first term on the RHS of Eq. (4.4) is just a precession, which 
keeps w + perpendicular to P + ; if we work with parallel transported axes we can ignore it. In 
the second term, the factor oi a/a means that the variation is significant only over a Hubble 
time. We then replace P + • P_ with its ensemble-averaged value, which is simply 2v 2 — 1. 
Again, this is not a fully controlled approximation, but should become progressively better 
for shorter segments. 

Along a left-moving characteristic we then have w + oc a 2v _1 oc 

r (2^-i)_ Here 

we take 

a = t r , so that r oc £ 1-r . Averaging over a translation-invariant ensemble of segments gives 

<[w + (a, t) - w + (a', r)] 2 > = r^-^)f{o - a') . (4.5) 
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The physical length I of a segment of coordinate length 5a is 

I = a(t)e(t)5a ~ a^f'^Sa , (4.6) 



where we have averaged in the e equation of motion ( 2.4[ )] | As time goes on, the physical 
length of the segment increases as a}~ 2v , but this is much slower than the growth of the 
Hubble length 0(t), so the segment 'propagates' to shorter and shorter scales. Correspond- 
ingly, as we evolve back, there will be a point at which its length approaches t and our 
approximations break down. At this point, we have to let the simulations deal with the 
horizon-scale dynamics, providing an initial condition that determine the constant t. That 



is, the correlator (4.5) takes some value when a 1 ~ 2v2 5a = Ct, where the constant C defines 
the matching scale; i.e. 5a = Ct 1+r ^ 2v2 ~" 1 \ This constant will be independent of time if the 
Hubble-length dynamics scales, so we can conclude that 

/(„-✓) = M |„ -of, X («) 



to cancel the time-dependence of the RHS of Eq. (4.5) at the matching point. In all, 



W i O", T — W_ 



.(a' 7 r)] 2 ) = 2A(l/t) 2x - (4-* 



Thus we obtain a definite form for the two-point function. Only a constant A need be 
taken from the simulations, aside from the ubiquitous v 2 . In particular, the latter determines 
the exponents Xm = 0.25 and Xr = 0.10. We compare the calculated two-point function with 
simulations |20j in Fig. 2. The log-log slope agrees rather well with the predicted value for 
a range of scales around near t, with A around 0.60 in both eras. At longer distance it goes 
over rather abruptly to zero slope, corresponding to the expected random walk. At shorter 
distances the simulations and the model differ. This is not a surprise. Our picture is that the 
structure propagates from the Hubble scale down to smaller scales. The expansion factor for 
the simulations [20] is only of order 3, so only a half-order of magnitude could have reached 
its asymptotic form; at shorter distances it must still reflect in part the initial conditions]^] 
In fact, the agreement is better than might be expected] 1 "] After a very long period of time, 
we would expect the power law to extend to arbitrarily short scales. Our exponents may 
be somewhat off, but the asymptotic power law form seems likely to be robust until new 
physics, particularly gravitational radiation, enters. 



The correlator of the tangent vector can be written in terms of Eq. (4.8), and the fact 



that this goes to zero with I means that the string become straighter at short distance. In 



To be precise, from the definition (2.3 1 of e it follows that I is the energy of the segment divided by fx. 
This is the length at an instant when the string is at rest, but it is convenient to work with this energy-length 
because it is conserved as the string oscillates. 

9 We note that the authors of Ref. [20] report that the same two-point function is obtained with other 
initial conditions, but this has not been quantified: it would be interesting in particular to compare the 
time-dependence with our form (14.51. 

ln There is an adage (Fermi?) that when one gets the physics right, one's approximations work better than 
they should. 
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X = 0.10 , A=0.60 




log W r q) 

Figure 2: Log-log comparison of the model (dashed line) with the data provided by Ref. 
in the radiation-dominated era (upper plot) and the matter-dominated era (lower plot). The 
x-axis is plotted in terms of a correlation length £, which is 0.30t in the radiation era and 
0.69t in the matter era. 
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other words, the fractal dimension goes to as 1 + 0([//t] 2 *). Thus we can answer the 
question that arose in the lensing discussion: long strings produce nice double images, with 
corrections [7] of order 1%. More complex images would be produced if we were near a (rare) 
large kink, or for lensing by a small enough loop. 

Most attempts to study cosmic string networks analytically have worked with a few 
macroscopic parameters, along the lines of the one-scale model. Ref. [6] is the most system- 
atic attempt to study the evolution of an ensemble analytically, but they did not attempt to 
explore the short- distance physics in a systematic way. Most of that paper is in the context 
of an approximation (their Eq. 2.24) that essentially assumes that the fractal dimension is 
larger than one. The small scale structure was treated by an Ansatz (their 3.16) which is 
equivalent to assuming that x — \- 



4.2 Very Small Scales 

At very small scales, we must take account of gravitational radiation, and also the matter- 



radiation transition. According to Eq. (4.6), the physical length of a segment grows as 



t r{\-2v ) ? Qr ^0.20 j n the matter 

era. Thus, a segment for which l/t ~ 1 at the time of matter- 
radiation equality will at later times have length l c {t) ~ t(t/t cq ) r( - 1 ~ 2v2 ^~ 1 ~ £(£/£ C q)~ ' 80 



Segments larger than this will have emerged from the Hubble scale during the matter era, 
and so their small-scale structure will simply be of the matter-dominated form. Segments 
smaller than this will have emerged during the radiation era, and so their structure will be 
of the radiation-dominated form with additional stretching during the matter era [TJ: 

<[w + (a) - w + ( ff ')f) = An(/A) 2Xm , / > k(t) , (4.9) 

= A(/ c (t)A) 2(xm ~ Xr) (/A) 2Xr = A c (t)(l/tf Xl , I < lc(t) . (4.10) 

In particular, the present value of A c is approximately 0.03. 

Now we can put in the smoothing via gravitational radiation, as discussed in Sec. 2.4. 
The details details of the derivation [TFJ [18j [8] are too lengthy for the present discussion, but 
we can summarize the result in a simple way. Denote the Fourier transform of the two-point 



function (4.8 ,4.10) by {w 2 (k)). Then, up to numerical constants, the string becomes smooth 
below a scale 

Z GW ~ TGfit(w 2 (k b )) , k b = 1/Gfit . (4.11) 
For the relevant tensions, the crossover form ( |4.10[ ) is the relevant one, and 

l GW ^TA c (t)(Gfi) 1+2x H . (4.12) 



11 This was also to conclusion of Ref. [20 ; the log- log plot accentuates the differences between our results 
and theirs, whereas the graphs of the fractal dimensions look more similar. 
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4.3 Loop Production 



We now take account of the production of small loops. Since this takes place on scales 
small compared to the Hubble scale, we can go to locally flat null coordinates u, v, where 
p+(u) = <9 u x and P-(v) = SyX, with u + v = t. 

A loop forms whenever a long string passes through itself, meaning that x(u, v + I) — 
x(u + l,v) for two points on the string. Defining 

L + (w,/) = / dup + (u) , L_(u, i)= / dvp-(v) , (4-13) 

J U J V 

this means that L+(u, I) = L_(t>,/) for some u, v, and I. The rate of loop formation, per 
unit volume in u,v,l, is then 

(det^tMx.O-MM))), J= V ^-^ M) . (4-14) 

Now, the components of L± are each proportional to the length I, so the 5-function implies 
a factor of l~ 3 . The rate would scale as Z~ 3 if the correlator of the w's were scale invariant. 



To work more carefully, separate p± as in Eq. (4.2), with the unit vectors P± proportional 
to L±. Then for small w±, 

i 3 (L + (ti,l) - i)) = rV(P + - P_)i(L + (M, I) - L(«, 1)) . (4.15) 

The magnitudes L± are both equal to 1(1 — 0(l 2x )), so the whole is of order Z~ 3-2x . The 
columns of J are of order l x ,l x ,l 2x (the last from the components parallel to P±), so the 
determinant is of order Z 4x giving l~ 3+2x for the rate. 

The total rate of string length going into loops is weighted by an additional factor of /, 



dir 2+2x . (4.16) 

This diverges at the lower end for x < ^, a crucial result [7]. This large production of small 
loops seems surprising at first because the string is becoming smooth in the sense that the 



fractal dimension approaches 1. What we have found in Eq. (4.16) is that the production 
of small loops is controlled by the rate at which the fractal dimension approaches 1, and in 
both the matter and radiation eras the approach is slower than the critical value \ = \- 



The total rate (4.16 ) is determined by energy conservation, Eq. (3.3 ). Naively this requires 



a lower cutoff on the integral, around I ~ O.lt. As we follow the evolution of the resulting 



loops, the same calculation leading to Eq. (4.16) shows that the production of small loops 



continues to diverge. This suggests a complicated fragmentation process. However, we have 
found a more physical way to introduce the cutoff [10J. If we separate the functions p± 
into a 'classical' long-distance piece and a random short distance part (with our calculated 
two-point function), we find that loop production takes place only near intersections of the 
long-distance parts of p + and p_ — that is, near cusps (Fig. 3). 
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We thank Carlos Martins and Paul Shellard for suggesting this possibility. 
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Figure 3: The functions p± written as a fixed long-distance piece plus a random short- 
distance part. Production of small loops takes place near the cusp of the long-distance 
part. 



Note that the figure is showing what is essentially the tangent vector to the string, rather 
than the string itself. The latter has a fractal dimension that approaches 1 at short distance, 
but the tangent vector path has a large fractal dimension as indicated. The two-point 
function (4.8) implies that a segment of length 5a makes a step of order Sa x . After N steps, 
the length of the path is N5a x , while the typical net distance covered (4.8) is (N5a) x , so the 
fractal dimension is 1/x- I n particular, this implies that near each cusp of the long-distance 
curve, their are very many small cusps. Formation of small loops occurs near these cusps: 
for a small loop, equality of L + and L_ occurs when the vectors p± are nearly parallel. 

Away from the cusp, production of loops is shut down, but as the cusp is approached 
we find that production of loops starts uniformly on all scales down to the gravitational 
radiation cutoff, rather than fragmentating from large loops down to small. This is still 
work in progress, but we appear to get a power law that extends all the way down to the 
gravitational radiation scale. 

It should also be noted that these small loops are moving rather rapidly [7] . The velocity 
is just L + /l, which by the above discussion is of order 1 — 0(A[l/t] 2x ), so the relativistic 7 
factor is of order A' 1 ^ 2 (l/t)~ x . Recalling from footnote [8] that / is defined in terms of the 
energy of a string loop, the rest mass of a loop is of order nlA l ^ 2 {l/t) x . 



19 



5 Conclusions 



5.1 A Synthesis? 

Some of the most recent simulations [SJ EJJ [22] have used a new technique to reach larger 
expansion factors than previously possible. At intervals, the configuration is simply copied 
eight times, to produce a periodic box of twice the size. A randomness is then briefly 
introduced into the evolution, to eliminate the unwanted long-distance correlations. Various 
tests appear to indicate that no significant artifacts are introduced by this process. 

These simulations show two peaks in the loop production. One is at a scale of order O.li, 
while the other remains near the short-distance cutoff. In Refs. j2Tj [22] it is argued that the 
smaller peak is an artifact of the initial conditions and will eventually disappear. However, 
it is this peak that we have found in the calculation of the previous section: these small 
loops are an inevitable result of the small scale structure on the long strings. We will be 
able to check this interpretation further when we complete the calculation of the small loop 
production function as described above. The distribution functions shown in Ref. [22] for 



the matter and radiation eras appear to match our Eq. (4.16) quite well, but a more careful 
study of loop fragmentation is needed [TO] . 

We have not found the large loops, but it is inevitable that some of these will arise from 
chance self- intersections of long strings. Our picture does explain why these large loops 
can survive fragmentation: one would expect from Fig. 3 that the regions near cusps would 
fragment extensively, but non-cusp regions of the long loop will remain. We must, at least 
for now, leave it to the simulations to determine the relative amount of string going into 
large and small loops. From Ref. [22] it appears that only 20% of the string goes into large 
loops during the matter era, and only 10% during the radiation era. However this reference 
also finds that the peak at large sizes is still growing at the end of the run, so the final 
numbers may be different. 

Thus it appears that we may be approaching an accord on the question of loop sizes. 
Some fraction of the string, perhaps a small fraction, goes into the largest scale under 
discussion, near the horizon. The rest goes into loops at the gravitational radiation scale, 
l{ ~ A c (t)(Gfi) 1+2Xl t. The simulations played an essential role in first revealing this double- 
peaked distribution, while the model plays a key role in confirming that the small loops, and 
in showing that gravitational radiation sets their actual size. 



The divergence (4.16) explains why there is not just a single scale in the string network. 
Indeed, the nonlinearities of string evolution have a remarkable ability to transfer energy 
from cosmic scales down to much shorter wavelengths. In this respect the cusps are like 
shock waves. Indeed, in shock waves the low energy field equations do break down, and 
shorter scale physics enters. This has also been suggested for string networks [251 [2EJ 127] - 
We believe that smoothing due to gravitational radiation prevents this, but cosmic string 
networks have been a source of many surprises, and perhaps more are in store. 
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5.2 Gravitational Wave Signatures and the Inverse Problem 



Now let us consider the signatures, assuming that the picture reached above is correct: most 
of the string length goes into loops near the gravitational radiation scale, and a fraction 



(10% in the radiation era and 20% in the matter era) goes into loops of size O.li. At pulsar 



frequencies, only the large loops are visible, so the energy density calculated in Sec. 3.4 is 



reduced by a factor of 10. The bound (3.7) is therefore weakened by a factor of 100, to 



G/i<2x 10" 7 , which is right around the CMB bounds (we have used the fact that these 
large loops are produced with relatively small velocities). 

In quoting such a "bound" it is important to note that the theoretical uncertainty in the 
network properties is not included, because this is not known. As one example, simulations 
of the full Abelian Higgs model seem to give long string densities that are a factor of three 
smaller [27] than simulations of the Nambu action, for reasons that are not understood. 
Including this uncertainty weakens the CMB bounds by a factor of 3, but the pulsar bounds 
by a factor of 9 because f2 G w sca l es as {Gijl) 1 / 2 . The pulsar bounds have a further uncertainty 
coming from the sizes and numbers of the large loops. 

The pulsar bound on Oqw is expected improve by an order of magnitude from further 
observations at the Parkes Pulsar Timing Array [39], improving the sensitivity to Gfj, by 
two orders of magnitude. The future Square Kilometer Array is expected to reach several 
orders of magnitude further in Oqw, but at some point a stochastic background from coa- 
lescing black holes will dominate the string signal. For the interferometers, Advanced LIGO 
sensitivity of dfl GW /d In z/ ~ 10 -8 to 10~ 9 implies sensitivity to G\i down to 10 -8 to 10~ 10 . 

Turning to the high harmonics of the large loops, Fig. 3 makes it clear that the large 
loops form without large cusps: the world-sheet region near the intersection of the classical 
p-t curves fragments into small loops. The functions p± have large discontinuities where 
the cusps would be. However, as the loops begin to decay via gravitational radiation, these 
curves smooth, because the higher harmonics decay faster than the lower. Thus the large 
loops will develop cusps as they decay. However, a factor of 10 is lost to the small loops 
(since large loops decaying today formed during the radiation era) so these are unlikely to 
be seen at LIGO except in very favorable parts of parameter space (such as small P). 

The small loops form at the scale set by the gravitational smoothing. We will adopt 
here our model from Sec. 4 to obtain representative numbers, but one should be mindful 
of the fact that the exponents and scales that we have found may be corrected by future 
work; we do not have a systematic way to estimate errors. For the relevant parameter range 



we will only see loops that form after the radiation era, so the relevant scale is Eq. (4.12) 



Zf ~ TA c (t)(Gfi) 1+2Xr t. However, because the loops are highly relativistic, their signatures 

13 Ref. 36J gives a more complete treatment of future observations, including possible backgrounds. Note 
however that it is based on a scenario where all of the string goes into loops of size O.lt, and so the actual 
signatures will be reduced by a factor of 10 according to our picture. Ref. ^3?j discusses gravitional waves 
from both large and small loops under varyious network scenarios. 
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are modified |8J. Thus, the formation scale l{ is \i 1 times the typical energy of a loop, giving 

eenergy-length ~ A c {t){G fif Xt . (5.1) 

This sets the period of oscillation as measured in the global FRW time, 

However, most of the emitted radiation is Doppler-shifted by a factor (1 — ~ ^y 2 ~ 
A c {t)~ l {I / t)~ 2Xl . For stochastic radiation the effective value of e is set by the period of this 
radiation, and so 

estoch ~ A(t) 1+2Xr r 2 *'(G>) 4 ^ 1+ ^ ~ 0.03(G» - 44 . (5.2) 



The frequency bound (3.10) becomes 



(GW L44 >-^- . (5.3) 

In the pulsar range this requires G\i > 3 x 10~ 7 , so these small loops will not be seen, but in 
the LISA range it is sufficient that G\i > 10 -11 , and in the LIGO range G\i > 2 x 10~ 14 . The 



density (3.13) becomes 

^ = 10c(^ )- 1/3 (G/i) - 52 . (5.4) 

In the LIGO band this is of order 3 x 10~ 6 c(Gyu) ' 52 , and in the LISA band it is of order 
10~ 4 c(G/z) a52 . Advanced LIGO sensitivity of dfi G w/^l nz/ o ~ 10~ 8 to 10~ 9 implies sensitiv- 
ity to Gfi only down to 10 -5 to 10~ 7 . Potential LISA sensitivity down to dQcw/dlnvo ~ 
10~ 12 would allow detection of tensions down to 10 -16 , but the frequency limit of 10~ n is 
encountered first. 

For the signal of individual cusps, on the other hand, the effect of the boost goes in the 
opposite direction: the effective size of the cusp is increased by a factor of 7, 



^cusp size ^ 



A(t) 1/2 ~ x T- 2 ^(G/i) Xr - 2 * ? . (5.5) 



This is around 0.02 over the range of interest. The effect of e cusp ~ 0.02 is to increase by a 
factor 0(5) the event rate for G\i greater than 10~ 9 , and decrease it by a factor 0(2) below 
this tension [42] . However, the relativistic loop velocities lead to two additional effects. First, 
the number of cusp events is controlled by e ene rgy- length, which is smaller by a factor of order 
50; the rate is increased by this factor. Second, the observed periodicity of a loop is set by 
e s toch : for G/i = 10~ 7 the periodicity is one year, but for Gfi = 10~ 8 it drops to two weeks. 
A one-year run will therefore see 25 cusps if it sees one, and so the mean event rate is not 



a good measure of the likelihood of seeing at least one event. Superimposing these effects 



on the results of Refs. jHJ H2J H3] suggests that these small-loop cusps are within reach of 



14 Thcre are other such statistical effects. For example, as Fig. 3 suggests, production of small loops is 
likely to be highly correlated in spacetime. Ultimately such effects may lead to improved sensitivity, by 
fitting to templates based on rare events. Ref. [S] examines one example of this type, the search for periodic 
signals from very small loops. 
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Advanced LIGO only for tensions very close to the current upper limit, but again it is clear 
that a more precise understanding of the networks is essential. 

Thus, of the four cases (large loops versus small and low harmonics versus high), the low 
harmonic stochastic radiation from large loops seems to be the most promising signature, 
with sensitivity to Gfi improving by several orders of magnitude in the next decade, both in 
pulsar observations and at Advanced LIGO. Assuming a vanilla string network, observation 
of the stochastic radiation would determine the string tension via (3.6), assuming that the 
large loop size and velocity distribution is accurately known, in particular the mean values 
of a 1 / 2 7~ 3 / 2 and the fraction of string going into long loops. 

In the approximation used in Sec. 3.4.1, the gravitational wave spectrum i/orfficw/^o is 
nearly flat. However, for Gfj, < 10 -8,5 , the pulsar observations begin to see loops decaying 
during the matter era (see the discussion below Eq. (3.5)), and so the energy density is 
larger, until G^i ~ 1CT 11 when the decaying loops are so small that their frequency is above 
the pulsar range. The slope of the spectrum (and, ideally, the observation of the stochastic 
radiation in both the pulsar band and the LIGO band) give some check on the cosmic string 
interpretation. 

It is interesting to go beyond the vanilla model and allow for a general value of P. The 
simplest argument would predict that the signals scale as l/P, because one would need to 
increase the long string density by this factor in order to get the same intercommutation 
rate. However, inserting a factor of P into the loop production in the one-scale model (2.7) 
one finds that the model predicts that the density scales as l/P 2 . This is because one needs 
to increase the rate of long string intercommutations, to produces more kinks, so that the 
loop production (which only comes from self-intersections) stays at the P = 1 value. On the 
other hand, simulations jUJHS] show a dependence more like P -0 - 6 ; evidently, the reduction 
of P is partly offset by multiple crossings of the colliding strings. 

There is some degeneracy in the determination of G\i and P. Note that the low frequency 
stochastic radiation ( 3.6[ ) and the high frequency stochastic radiation (5.4) depend on Gfi 
approximately as (G/i) 0,5 , and so on Gfi and P only through (Gfi)°' 5 P~ 0S . Thus, measure- 
ment of the overall normalization of the stochastic radiation does not break the degeneracy 
between P and Gfi even if both the low harmonic and high harmonic contributions are seen 
(unless Gfi < 10~ 8 ' 5 , as discussed above). 

Ultimately, LISA will be sensitive to low-harmonic gravitational radiation from large 
loops, and individual cusps from large and small loops, over the most or all of the range of 
tensions that arise in brane inflation. This will allow several independent determinations of 
Gfi and P, making it possible to test the cosmic string hypothesis, and look for evidence of 
a less vanilla type of cosmic string. 
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